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ABSTRACT 

We compare estimates of stellar mass, M*, and dynamical mass, Ma, for a sample of galaxies from the Sloan 
Digital Sky Survey (SDSS). Under the assumption of dynamical homology (ie., ~ v^Re, where cto is the 
central velocity dispersion and R e is the effective radius), we find a tight but strongly non-linear relation be- 
tween the two mass estimates: the best fit relation is M* oc M° 73 , with an observed scatter of 0.15 dex. We also 

find that, at fixed M*, the ratio M*/Md depends strongly on galaxy structure, as parameterized by Sersic index, 
n. The size of the differential effect is on the order of 0.6 dex across 2 < n < 10. The apparent n-dependence of 
M*/M,j is qualitatively and quantitatively similar to expectations from simple, spherical and isotropic dynam- 
ical models, indicating that assuming homology gives the wrong dynamical mass. To explore this possibility, 
we have also derived dynamical mass estimates that explicitly account for differences in galaxies' structures. 
Using this 'structure-corrected' dynamical mass estimator, Mj „, the best fit relation is M + oc ^ ) ^ 2±0 01 (±o o8) 
with an observed scatter of 0.13 dex. While the data are thus consistent with a linear relation, they do prefer 
a slightly shallower slope. Further, we see only a small residual trend in M*/Md,„ with n. We find no statisti- 
cally significant systematic trends in M*/M± n as a function of observed quantities (e.g., apparent magnitude, 
redshift), or as a function of tracers of stellar populations (e.g., Ha equivalent width, mean stellar age), nor 
do we find significantly different behavior for different kinds of galaxies (ie., central versus satellite galaxies, 
emission versus non-emission galaxies). At 99 % confidence, the net differential bias in M*/Md,« across a wide 
range of stellar populations and star formation activities is < 0.12 dex (« 40 %). The very good agreement 
between stellar mass and structure-corrected dynamical mass strongly suggests, but does not unambiguously 
prove, that: 1.) galaxy non-homology has a major impact on dynamical mass estimates, and 2.) there are not 
strong systematic biases in the stellar mass-to-light ratios derived from broadband optical SEDs. Further, ac- 
cepting the validity of both our stellar- and dynamical-mass estimates, these results suggest that that the central 
dark-to-luminous mass ratio has a relatively weak mass dependence, but a very small scatter at fixed mass. 
Subject headings: galaxies: fundamental parameters — galaxies: kinematics and dynamics — galaxies: stellar 
content — galaxies: structure 



1. INTRODUCTION 

Techniques for estimating galaxies' stellar masses are a 
crucial tool for understanding galaxies and their evolution. 
There are tight and well-defined correlations between stellar 
mass and many other important global properties like color, 
size, structure, me tallicity, star formation activity, and envi- 
ronment (see, e.g., [K auffmann et al. 2003b; Shen et alj|2003b 
iBlanton etaL 2005b: iGallazzi et ^1.1120061) 7 Given a galaxy's 
stellar mass, M*, it is thus possible to predict a wide vari- 
ety of global properties with considerable accuracy. In this 
sense, stellar mass appears to be a key parameter in determin- 
ing (or at least describing) a galaxy's current state of evolu- 
tion. Moreover, since the growth of stellar mass (cf. absolute 
luminosity, color, etc.) is relatively slow and approximately 
monotonic, stellar mass is a particularly useful parameter for 
quantifying galaxy evolution. 

Stellar mass estimates, whether derived from spectro- 
scopic or photometric spectral energy distributions (SEDs), 
are plagued by a variety of random and systematic errors. 
These include a generic degeneracy between mean stellar age, 
metallicity, and dust obscuration. It is typical to make the sim- 
plifying assumptions that galaxies' stellar populations can be 
described en masse (ie., neglecting age/metallicity gradients 



and complex dust geometries), and that galaxies' complex star 
formation histories can be described parametrically. It is rare 
to attempt to account for active galactic nucleus (AGN) emis- 
sion. The stellar initial mass function (IMF), including its 
universality or otherwise, remains a major 'known unknown'. 
Then there is the complication that different wavelengths 
probe different aspects of the stellar population; the inclusion 
of restframe UV or NIR data can thus, in principle and in 
practice, have a large impact on the estimated stellar mass. 
These effects are compounded by uncertainties in the stellar 
evolution models themselves. A topical example is the impor- 
tance of NIR-luminous thermally pulsating asymptotic giant 
branch (TP-AGB) stars: for t he same data and stellar p op- 
ulation paramete rs, the use of Bruzual & Chariot] (1200 31) or 
Maraston (2005) models can change the derived value of M* 
by a factor of 3 for galaxies that host young (< 1 Gyr) stars, 
but on ly if restframe NIR data are included ( van der Wel et al. | 
120061; iKannappan & Gawiserl 120071) . iConrov. Gunn & White! 
(2009) have argued that the total random uncertainties in M* 
are on the order of ~ 0.3 dex for galaxies at z ~ 0. 

For these reasons, it is essential to devise some way of 
assessing the quality of stellar mass estimates through com- 
parison to some other fiducial mass estimate — this is the 
primary motivation for the present paper. Specifically, us- 



2 



Taylor, Franx, Brinchmann, van der Wei, van Dokkum 



ing a number of the publicly available 'v alue added' cata- 
logs of the Sloan Digital Sky Survey (SDSS: lYork et al.ll2000t 
Straus s~et alj2002l ). we will compare stellar mass estimates to 
total mass estimates derived from galaxy dynamics. 

From the outset, we note that a difference between two 
quantities shows only that: a difference. With no definitive 
standard to use as a basis for comparison, the best that we 
can hope for is consistency between the two mass estimates. 
Further, if and when there are differences, it is impossible to 
unambiguously identify where the 'fault' lies — or even if 
there is indeed a fault. For example, it is likely that the ratio 
between stellar and total mass varies as a function of mass, 
and/or some other global property/ies. We will also explore 
this issue in some detail. 

This kind of comparison has been done for SDSS galax- 
ies by iDrory. Bender & H opp (2004), who considered both 
stellar mas s estimates derived from t he SDSS spectra, as de- 
scribed by Kauffmann et al. (2003a), and those derived us- 
ing SED-fitting techniques that are commonly used at for 
higher-redshift studies. These authors find a relatively tight 
correlation between the two stellar mass estimates, with 
a mild systematic bias depending on Ha equivalent width 
(EW). This bias suggests a potential problem with the stel- 
lar mass estimates as a function of specific star forma- 
tion rate (SSFR). Further, both stellar mass estimates cor- 
related well with the simple dynamical mass estimate, Md 
(defined below), but showed a clear trend in M*/Md with 
mass, such that less massive galaxies had higher values 
of M*/M d . IPadmanabh an et all ([2004) used the mass de- 
pendence of M»/Md to argue for an increasing stellar-to- 
dark mass ratio for ellipt ical gala xies with higher m asses, 
as did iGallazzi et all (12005b . Both iRettura et all (12006b and 
Ivan der Wei et al.l (1200 6) have performed similar comparisons 
for z < 1 galaxies, with similar conclusions. 

It is common practice to derive a simple dynamical mass 
estimate based on the velocity dispersion, <to, and the effective 
radius, R & , via the scalar virial theorem: 



' k<7%R e 



(1) 



(This is the dynamical mass estimator used for each of the 
studies cited in the previous paragraph.) The constant k is 
usually assumed to be in the range 3-5, and is intended to ac- 
count for the 'degree of virialization', including the effects of 
dark matter and t he intrinsic shape of the velocity dispersion 
profile (see. e.g..)Cappellarietal.l l2006; Gal lazzi et al.ll2006t 
Ivan der Wei et al.l 120061) . By assuming a constant k for all 
galaxies, this expression implicitly assumes that all galaxies 
are dynamically homologous, or self-similarQ 

But it is important to remember that the observed veloc- 
ity dispersion is actually the luminosity-weighted mean of the 
true, radially-dependent velocity dispersion, projected onto 
the line of sight, and within the spectroscopic aperture. The 
shape of mass profile has a strong influence on the spatial and 
dynamical distribution of stellar orbits: in general, the relation 
between the observed velocity dispersion and the underlying 
mass profile thus depends on structure as well as size. As a 
dynamical mass estimator, M& can therefore only be consid- 
ered as approximate. (The tilde in M& is intended to remind 



1 But see also Wolf et al. 12009), who derive a mathematically identical 
relation from the spherical Jeans equation for a system in dynamical equilib- 



th e reader of this fact.) 

iBertin. Ciotti & Del Principel ( l2002l) provide an analytic ex- 
pression that makes it possible to approximately account for 
this effect. Using their formulation of the problem, the dy- 
namical mass can be expressed as: 



GM A . n = K v (n)at ) R e 



(2) 



Here, the term Ky{n) encapsulates the effects of structure 
on stellar dynamics. (The subscript n in Mj „ is intended to 
make it clear that non-homology has been accounted for as a 
function of Sersic index, n.) For convenience, we will refer to 
Md „ as a 'structure corrected' dynamical mass estimator, but 
we note that the inclusio n of a structure-dependent t erm is 
not strictly a correction. Bertin, Ciotti & Del Principd d2002l) 
also provide an analytical approximation for Kv(n): 



K v (n) S 



73.32 



10.465 + («-0.95) 2 



: + 0.954 



(3) 



This expression for Kv(n) has been derived assuming a 
spherical mass distribution that is dynamically isotropic and 
non-r otating, and which, in projection, follows a Sersic] dl963l 
1968) surface density profile. For this (admittedly simple) 
scenario, this approximate expression for Ky{n) is accurate at 
the percent level for 1 < n < 10. Substituting trial values of 
n = 2 and n = 8 into Equation [3] suggests that the differential 
effect of non-homology on the inferred value of the dynamical 
mass is as much as a factor of 3, or 0.5 dex. Our first task in 
this paper, then, will be to explore the importance of structural 
differences between galaxies, using this prescription. 

Before we begin, note that there are alternative approaches 
to exploring the consistency between stellar and dynamical 
mass estimates. In particular, a number of authors have con- 
sidered the relation between galaxies' stellar and dynamical 
masses in the context of well known scaling re lations between 
lumino sity /mass and dynamics. For example. iBell & de Jongl 
d2001l) considered the relation between baryonic (cf. stellar) 
mass, Mbar, and circular rotation velocity, Vc, for disk galaxies 

— the baryonic Tully-Fisher relation. These authors showed 
that stellar mass estimates based on different passbands (ie. 
M*/Lv versus M*/Lk) produced consistent Mbar— Vc relations. 
Furthermore, for a fixed IMF, they argued that it was possible 
use a single color to estimate stellar mass-to-light ratios with 
an accuracy of 0.1-0.2 dex. 

There have also been a number of analogous studies 
for elliptical galaxies, based on the fundam ental plane 
(Diorgovs kv & Davislll987l : iDressler et alj [l987). which can 
be understood as a correlation between the dynamical mass- 
to-light ratio, M&/L, and surface brightness. These stud- 
ies (see, e.g., Cappellari et al. 2006; and references therein) 
have tended to focus on the 'tilt' of the fundamental plane 

— that is, the deviation of the observed relation from the 
expectation assuming both a constant M*/L and structural 
and dynamical homology for all early type galaxies. The 
tilt of the fundamental plane thus offers a means of prob- 
ing variations in M&/L (including both variations in M»/L 
due to different stellar populations, and variations in M*/Md 
due to, e.g., different dark-to-luminous mass ratios) and/or 
the degree of non/homology. While the relative impor- 
tance of these different effects remains an open question, 
it seems clear that, at least for early type galaxies, both 
Af»/L and M^/Md vary systematically with mass (see, e.g. , 
Prugn iel & Simienll 19961: IBertin. Ciotti & Del Principell2002l: 
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Truiillo et al. 2004; 


Cappellari et al.ll2006l lLa Barbera 2008; 


Allanson et al. 2009 


). 



This paper is structured as follows: in Section [2] we 
describe the various SDSS-derived catalogs that we will use, 
including the definition of our galaxy sample. We validate the 
velocity dispersion measurements used to derive dynamical 
masses in Appendix [A] In Sections [3] and H] we present 
two parallel comparisons between stellar and dynamical 
mass estimates for galaxies in our sample. First, in Section 
[3] we use the simple dynamical mass estimate then, 
in Section [4] we show how the comparison changes using 
the structure-corrected dynamical mass estimate, M& M . In 
Section [5] we explore the consistency between M&„ and 
M*. In particular, we will show that there are no statistically 
significant trends in the ratio M*/Md,„ that would indicate 
measurement biases in M» and/or M&„; this is not the case 
for the simple dynamical mass estimate We show in 
Appendix [B] that these results are not unique to the sample 
we consider in the main text. We discuss the interpretation 
and implications of this result in Section [6] before providing 
a summary of our main results and conclusions in Section|7] 

Throughout this work, we will assume the concordance cos- 
mology; viz.., (O m , » a, n Q ) = (0.3, .7, 1.0), and H = 70 
km/s/Mpc, and adopt a lChabrierl (120031) IMF. 

2. DATA 

This work is based on data drawn from several publicly 
available catalogs based on the SDSS dataset. Our anal- 
ysis is based on redshifts and ve locity dispersions fro m 
the basic SDSS catalog for DR7 dAbazaiian et al.1 12 009H 
We u se Sersic-fit structural parameters from iGou et al] 
(120091 hereafter lG09h and SED-fit stellar mass-to-light 
ratio measurements from the DR7 Max-Planck-Institute for 
Astrophysics (MPA)/Johns Hopkins University (JHU) value 
added catalog In Appendix IB~1 we repeat our analysis using 
the Sersic -fit structural parameters given in the New York 
University (NYU) V alue Added Galaxy Catalog (VAGC; 
iBlanton et al.ll2005ah for DR7. Each of these catalogs have 
been well described and documented in the references given; 
in this Section, we only briefly summarize the most relevant 
aspects of each catalog for the present work. 

2.1. Redshifts and Velocity Dispersions 

There are two sets of redshift and velocity dispersion mea- 
surements given in the basic SDSS catalog for DR6 and DR7: 
the 'spectrolD' values produced by the Chicago group, and 
the 'specBS' values produced by the Princeton group. In 
terms of redshifts, the two algorithms produce virtually iden- 
tical results. The major difference between the two algorithms 
is that, whereas the Chicago pipeline only gives velocity dis- 
persion measurements to those galaxies that are spectrally 
classified as being 'early type', all galaxies are given a veloc- 
ity dispersion measurement by the Princeton pipeline. From 
DR6, both the Princeton and Chicago velocity dispersion 
measurement algorithms have been updated, so as to eliminate 
the sy stematic bias at low dispersions identified by iBernardil 
(2007) for the DR5 valuesQ In Appendix |A] we compare 

2 Accessed via t he Catalog Archive Server (CAS; IThakar et aTll2008l) : 

http://casjobs.sdss.org/CasJobs/] 

1 Available via http://www.mpa-garching.mpg.de/SDSS/ 
4 See http://www.sdss.org/dr7/algorithms/veldisp.html for a discussion of 
the spectrolD and specBS algorithms, as well as a comparison between these 
values and those from Bernardi et al. 1 2003a b) and SDSS DR5. 



both se ts of velocity dispersions to those given by Faber et al. 
( 1989) for brig ht, early typ e galax ies: in both cases, the values 
agree with the lFaber et al.l dl989l) catalog values with an rms 
difference of ~ 18 km/s and no discernible systematic bias. 

The default redshifts and velocity dispersions for in the 
SDSS catalog (specifically, using SDSS parlance, the pa- 
rameters z and veldisp given in the table specOb jAll) 
are the Chicago values. For the sole reason that Princeton 
velocity dispersions are given for all galaxies (rather than 
only the spectroscopically early-types) we have elected to use 
the Princeton values instead; these are also the measurements 
adopted for both the NYU and the MPA/JHU value added 
catalogs. 

As we have mentioned in the Introduction, the observed 
velocity dispersion is the luminosity weighted average within 
the (projected) spectroscopic aperture. In order to account 
for aperture effects, we have scaled the observed value, er b, 
for each galaxy to a central velocity dispersion, <to, which is 
defined to be that that would be observed within a circular 
aperture wi th a radius equal to 1/8 times the apparent e ffectiv e 
radius, 6 e Jcir gensen. Franx & Kjaergaardl (see, e.g. 119951) . 
This correction has been made assuming a(R) oc 066 ; i e _ 7 
o-o/CTob = (8e a p/e e )- ' 066 , where 9 ap = l."5 is the radius of 
the SDSS spectroscopic aperture. The scaling of a(R) has 
been derived by taking a luminosity-weighted integral of the 
spatially resolved velocity dispersions o f galaxies from the 
SAURON survey dCappellari et alj 120061) . The corrections 
themselves are small — the median correction is 0.02 dex, 
with an rms scatter of 0.02 dex — and does not have a major 
impact on our results. Our qualitative conclusions do not 
change if we assume the slightly weaker radial dependence 
a(R) o c R~° 04 as found by lJ0rgensen. Fran x & Kiaergaard 
(fT995h . or if we neglect this correction altogether. 

2.2. Sersic Parameters: Size, Flux, and Structure 

IG09I have derived r-band structural parameters including 
total magnitude, ;« tot , effective radius, O e , and Sersic index, «, 
for a modest sized sample of SDSS gala xies. (We discuss the 
specific sample selection in Section l2~4l below.) These values 
have been derived via parametric fits to the (2D) r-ba nd sur- 
face brightn ess distribution of each galaxy, assuming a lSersid 
d 1 9631 fl968b profile, and convolved with th e appropriate PSF , 
using the publicly available code galf it (Peng et al. 2002). 
In order to account for blending, where two galaxies are very 
close, both the target and companion(s) are fit simultane- 
ously. Through analysis of simulated data, IG09I show that the 
median error in each of m tot , O e , and n to be less than 10 %. 
(Although at the same time, they show that the uncertainties 
on the derived parameters associated with background 
subtraction alone can be significantly greater than this.) 

In Appendix iBl we will also make use of Sersic-fit struc - 
tural parameters from the NYU VAGC dBlanton et alJl2005d) . 
Whereas galf it considers the 2D surface brightness distri- 
bution, the VAGC algorithm makes fits to the ID azimuthally 
averaged cur ve of growth. The a nalysis of simulated data 
presented by Blant on et alJ ((2005a) shows the VAGC Sersic 
parameters to be systematically biased towards low fluxes, 
sizes, and Sersic indices. This problem becomes progres- 
sively worse for larger n, such that sizes are underestimated 
by > 20 % and fluxes by > 10 % for n > 5. lG09l have shown 
that this bias is produced by background over-estimation and 
over-subtraction in the VAGC Sersic fits, owing to the use of 
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a 'local', rather than a 'global' background estimator. 

2.3. Stellar Masses 

We note that there are rather larg e diff erences between 
the the Sersic magnitudes given by IG09I and the default 
model magnitudes given in the SDSS catalog. The model 
photometry is derived by making parametric fits to the 2D 
surface brightness distribution in each band, usin g the sector 
fitting technique described by Stra uss et al.l (|2002). These fits 
assume either an exponential or a De Vaucouleurs profile; the 
profile shape is chosen based on the best-fit \ 2 m the r-band. 
For galaxies that are best fit by a De Vaucouleurs model we 
find (»?G09,r _ '«DeV.r) ~ — 0.26 + 0. 1 l(n — 4), where n is the 
Sersic index reported by G09; the scatter around this relation 
is at the level of 0.15 mag (la). That is, even where G09 
find n = 4, their flux is approximately 0.26 mag brighter than 
the SDSS (De Vaucouleurs) model flux; this discrepancy 
is larger for larger values of n. Fo r this reason, we take the 
r-band Sersic magnitude from G09 as a measure of total flux. 

To derive a stellar mass, we then use M*/Ls taken from 
the MPA-JHU catalog (DR7), which is maintained by the 
Garching group|_ N ote that, unlike previous MPA-JHU cat- 
alogs (e.g. Kauffma nn et al.l l2003a; Brinch mann et al.ll2004t 
Galla zzi et al.l 120051) . which were based on the SDSS spec- 
troscopy, these masses are derived from fits to the ugriz 
model SEDsQ Note, however, that the SED photometry has 
been corrected for emission lines, according to the line-to- 
continuum flux ratio in the spectroscopic fiber apert ure (we 
discuss the importance of this correction in Section 15.31 be- 
low). The SED fits are b ased on the synthetic stellar popu- 
lation library described by Gallazz i et al.l (120051). whic h have 
been constructed using the iBruzual & Charlotl (120031) stellar 
population models, assuming a Chabrier (2003) stellar IMF. 
These M*/Ls have been shown to be in good agreement (rms 
in AlogM* of < 0.1 dex) with the sp ectrally-derived values 
described bv lKauffmann et al.l (I2003al) for DR4_] 

2.4. Sample Definition 

Our decision to rely on the Sersic structural parameters 
derived by G09 restricts us to working with their sample. Our 
rationale for this decision is that, in our estimation, the G09 
Sersic fits are the most robust that are presently available. 

The G09 sample was constructed with the specific goal 
of exploring differences in the sizes and structures of 'cen- 
tral' and 'satellite' galaxies in groups and clusters. To this 
end, they selected 911 z < 0.08 'centrals' as the first-ranked 
(in ter ms of M„) group/cluster members from the lYang et al.l 
( 120071) group catalog, which was in turn constructed from the 
DR4 NYU VAGC. These galaxies were selected to have a flat 
logarithmic distribution in halo mass in the range 11.85 < 
logM ha io/M Q < 13.85 (800 galaxies), plus 100 galaxies in 
the range 13.85 < logM halo / M < 14.35, and all 11 central 
galaxies in clusters with logMh a io/M0 > 14.35. In this way, 
the central galaxy sample was constructed to span a represen- 
tative range of (large) halo masses. 

G09 also construct two z < 0.08 'satellite' control samples, 
in which the satellite galaxies are selected to match the central 

5 Available via http://www.mpa-garching.mpg.de/SDSS/DR7/ 

6 Note that in the SDSS algorithm, when deriving the ugiz model photom- 
etry, the structural parameters in the fits are held fixed to the r-band values; 
only the overall normalization (ie., total flux) is allowed to vary. The fits 
in each band are also convolved with the appropriate PSF. In this sense, the 
model SEDs are both aperture matched, and PSF-corrected. 

7 See http://www.mpa-garching.mpg.de/SDSS/DR7/mass_comp.html 



galaxies. For the first of these, satellites are chosen to match 
centrals in M* to within 0.08 dex; in the second, satellites are 
also required to match centrals to within 0.03 mag in (g—i) 
color. Because more massive galaxies are more likely to be 
(counted as) centrals, not every central has a satellite counter- 
part within these limits: the matching is more than 90 % suc- 
cessful forM* < 10 10 85 M , and less than 10 % successful for 
> 10 1115 M Q . The two satellite samples, so constructed, 
consist of 769 and 746 galaxies, respectively. 

G09 exclude a number of these galaxies from their analysis 
because of confusion, leaving a sample of 879 central 
galaxies, and two samples of 704 and 696 satellites each. 
While duplicates are not allowed within the individual 
satellite samples, some galaxies do appear in both samples; 
combining the two satellite samples we have 1167 unique 
galaxies. We exclude a further 71 galaxies whose spectra 
are not deemed 'science worthy' by the SDSS team (ie. the 
flag sciencePrimary is set to zero). In order to avoid 
very large errors in the dynamical mass estimates, we also 
exclude 160 galaxies that have relative errors in their velocity 
dispersion measurements that are greater than 10 %. This 
requirement excludes mostly low-n and low-M* galaxies: the 
vast majority of these 160 galaxies have M* < 10 10 8 M Q and 
n < 1.5. Our results do not depend on these selections. We 
are thus left with a sample of 1816 galaxies, of which 784 
have been selected as central galaxies, and 1032 have been 
selected as being satellites of comparable mass. 

The major disadvantage to using the G09 sample is that 
the relative number of central/satellite galaxies is not at all 
representative of the general galaxy population. However 
it is worth noting that G09 have shown that, at least for 
structurally early type galaxies, and after matching both color 
and mass, there are no differences in the sizes and structures 
of central and satellite galaxies. They thus conclude that the 
distinction between central and satellite has no impact on 
galaxy structure. This already suggests that the G09 sample 
may be adequately statistically representative of the massive 
galaxy population. 

Even so, we will explicitly examine the possible role of 
sam ple selection effects in shaping our results in Section 
I5.1l bv comparing different subsamples from within the com- 
bined G09 sample. Further, in Appendix [B] we analyze 
a more general galaxy sample, using structural parameters 
from either the NYU VAGC Sersic fits or the SDSS De Vau- 
couleurs/exponential model fits. 



3. RESULTS I — COMPARING STELLAR AND 
DYNAMICAL MASS ESTIMATES ASSUMING 
DYNAMICAL HOMOLOGY 

In this Section and the next, we present parallel compar- 
isons between stellar mass and two different estimates of dy- 
namical mass. As we have said in the Introduction, it is 
common practice to obtain a simple dynamical mass estimate 
based on oo and R e alone, using the scalar virial theorem; viz. 
M& ss ka^R^. In Section |3~T1 we directly compare the values 

of M* and for the IG09I galax ies; we will assume k = 4. 
We will then argue in Section l3~2l that the agreement between 
stellar and dynamical mass estimates may be significantly im- 
proved if we allow for non-homology. To test this idea, in 
Section |4j we will perform the same comparisons using the 
structure corrected dynamical mass estimator, Mj,,,. 
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FIG. 1. — Comparing stellar and dynamical mass estimates under the as- 
sumption of dynamical homology. — This Figure shows the relation be- 
tween stellar mass and a simple estimate of dynamical mass, GMj = 4<Tgi? c , 
for galaxies in the G09 sample. The black points show the data themselves; 
the red points with error bars show the median and 16/84 percentile values 
of M* in narrow bins of Mj ; the heavy dashed line shows a log-linear fit to 
the data, with the form and parameters as given. While there is a relatively 
tight correspondence between M* and M&, the relation is clearly non-linear. 
Further, for < 10 10 ' 5 Mq, M* appears to exceed M L \ for the majority of 
galaxies, which is logically inconsistent. We explore these results in greater 
detail in Figures |2]and[3] In Figure E we show how these results change if 
we account for structural and dynamical non-homology in our estimates of 
dynamical mass. 

3.1. The Relation Between Stellar and Dynamical Mass 

' In Figured] we compare the values of the simple dynami- 
cal mass estimator, M&, to the values of M* for galaxies in the 
G09 sample. The first thing to notice is that there is a rela- 
tively tight but clearly nonlinear relation between M d and M*, 
such that M* oc M d with a < 1 . Moreover, this simple analy- 
sis suggests that for many galaxies, including the majority of 
galaxies withM d < 10 10 5 M Q , M* actually exceeds M A . This 
is logically inconsistent, and necessarily implies a problem in 
the calculation of M* and/or M d . 

Before discussing this result further, it is appropriate to 
make a few comments about the random errors in our esti- 
mates of M* and M d . In particular, it is important to realize 
that the errors in O e , m tot , and n are strongly covariant: for 
example, an error in the structural index will affect the values 
of both e and m tot . Because M» depends on the measured 
value of m tot , and M d on the measured value of 9 e , M* and 
Ma are thus also covariant. This makes the seemingly trivial 
task of fitting a line to the observed M„-M d relation rather 
problematic. To do this properly would require full and con- 
sistent treatment of the covariant uncertainties i n the S ersic-fit 
parameters, but this information is not given bv lG09L 

Our solution to this problem is simply to minimize the mean 
absolute perpendicular distance between the data and the fit. 
When doing so, we also use a 'sigma-clipping' algorithm to 
avoid the influence of the most egregious outliers; specifi- 
cally, we iteratively exclude points that lie off the best-fit line 
by more than 5 times the rms offset. While the gradient of 



the best-fit line does depend on the fitti ng sc heme used (we 
will explore this in more detail in Section l4Tt . the best fit pa- 
rameters are not strongly dependent on how aggressively we 
sigma-clip. In order to avoid strong covariances between the 
slope and intercept of the best-fit line, we actually compute 
the fit in terms of log(M d /lO n M0); that is, we fit a relation 
of the form y = a(x- 1 \) + b\\. Statistical uncertainties on the 
fit parameters have been derived from bootstrap resampling. 
The best fit to the M*-M d relation, so derived, is shown as 
line heavy dashed line in FigureQ] The best fit parameters are 
a = 0.73 ±0.007 and b n =-0.14 ±0.003. 

In Figure |2] we explore the relation between M* and M A 
in greater detail. The different panels of this Figure show 
the difference between M* and M d as a function of several 
interesting global properties: (from left to right) galaxy 
mass, structure, dynamics, and surface density. It clear that 
that M»/M d is strongly correlated with all four of these 
parameters. For each of the parameters shown, the size of 
the median trend in M»/M d across the sample is on the order 
of 0.5 dex, although it is slightly lower for M* and slightly 
higher for effective surface density. 

To quantify this statement, we have again made fits to the 
data, assuming the form y = a(x — X) + bx, where X is an 
arbitrary value chosen to be close to the median value of 
the quantity x for our sample. For these fits, in contrast to 
the previous Section, we have minimized the mean absolute 
vertical offset between the data and the fit. Again, we use 
a non-aggressive sigma-clipping scheme to exclude extreme 
outliers. (In all that follows, when considering the stellar- 
to-dynamical mass ratio, we will always fit in this way; we 
will only use the minimum perpendicular distance algorithm 
described above when fitting the relation between stellar and 
dynamical masses.) The best fit lines to the data, so derived, 
are shown as the heavy dashed lines in each panel; the best-fit 
parameters are given in each panel. In the case of Sersic 
index, the scatter around the best fit relation is « 0.12 dex; for 
the other three parameters it is slightly higher: rs 0.15 dex. 

3.2. The Importance of Galaxy Structure 
in Dynamical Mass Estimates 

There are of course strong correlations between mass, 
velocity dispersion, surface density, and structure. It is thus 
possible that the apparent trend with any given parameter 
in Figure |2] is 'spurious', in the sense that it is driven by a 
trend in another more 'fundamental' parameter. We note 
that galaxies' star formation activity and histories have 
been shown to correlate closely with all of mass, velocity 
dispers io n, and surface de nsity (see, e.g., [K auffman n et al. 
| 2003allR I2006t IFranx etaLl 120081; iGraves, Faber & Schiavon 
2009). Indeed, with the assumption that M* ~ M d .„, these 
three quantities are all related by factors of R e , which is also 
closely correlated with M* dShen et al.l 120031; IFranx et ail 
2008). But the fact that M»/M d depends on galaxy structure 
— and particularly the agreement between the observed trend 
and expectations derived from a simple dynamical model — 
immediately suggests that structure-dependent differences 
in galaxy dynamics may play a role in the results shown 
in Figure Q] With this in mind, in Figure we attempt 
to separate out the M„- and n-dependences of M*/M&. 
Specifically, we want to test the hypothesis that departure 
from linearity in the M»-M d relation seen in Figure Q] is at 
least in part a function of structure, and not mass. 
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Fig. 2. — Comparing stellar and dynamical mass estimates under the assumption of dynamical homology. — Each panel of this Figure plots the stellar-to- 
dynamical mass ratio, M*/Md, as a function of (left to right) mass, structure, velocity dispersion, and effective surface density. Within each panel, the black 
points show the data themselves; points that fall outside the plotted range are shown with a small grey plus; the large points with error bars show the median 
and 16/84 percentiles of log(M»/M,j) in narrow bins of the quantity shown on fhejc-axis. The dashed lines show fits to the data of the form y = a(x— X) + bx, in 
analogy to FigurefJJ At least when using this (overly) simple way of estimating galaxies' dynamical masses, there are strong trends in M t /M^ with both mass and 
structure. We see similarly tight and strong trends with velocity dispersion and effective surface density, as well as with other parameters like size (not shown). 
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FIG. 3. — Separating out the mass- and structure-dependences of the stellar-to-dynamical mass ratio, M»/Md. — In the left panel, we plot M*/Mi as a function 
of M f ; the thick lines in this panel show the median relation in bins of n. In the right panel, we do the reverse: M* /M& is plotted as a function of n, with the solid 
lines showing the median relation in bins of M* . In both panels, we give the parameters of the best fit relation for galaxies in each of the bins shown. It is clear 
from both panels that affixed M t , the scatter in M,/Ma is directly related to n. It is also true that, even at fixed n, M f /Ma appears to vary with M f ; we explore 
this issue further in Figure|6] For this Figure, we have assumed = 4o^R c ; in the right-hand panel, the solid curve shows the expected shape of the M*/M$ 
relation for a spherical and dynamically isotropic system that follows a Sersic profile, derived using Equation [5] The general agreement between the shape of 
this curve and that of the M*jM&-n relation suggests that including this term may significantly improve the correspondence between stellar and dynamical mass 
estimates. We explore this issue further in Figures 



Figure [3^ shows M^/M^ as a function of M*; the colored 
lines show the median relation in bins of Sersic index. The 
median relation between M*/Md and M* has a rather similar 
slope for each different n bin: M»/Md does depend on mass. 
If the dynamical mass-to-light ratio were to depend on mass 
only, however, we would expect the relations for different 
Sersic indices to overlap. Instead, the relations for each bin 
are clearly offset from one another. That is, at fixed mass, the 
scatter in M*/Md is closely correlated with galaxy structure. 

In Figure Qj), we do the opposite: in this panel, we plot 
M*/M& as a function of Sersic index, and the different lines 
show the median relation in bins of stellar mass. Again, it 
is clear that M*/Md depends on both M* and n: the median 
relations for each different mass bin are roughly parallel, but 



offset from one another. 

Further, the different mass bins in Figure|3}) cover different 
ranges in Sersic index: where the lowest mass bin is dom- 
inated by galaxies with 1 < n < 4, the highest mass bin is 
dominated by n > 4 galaxies. Similar behavior can be seen 
in Figure [3^: the lowest n bin contains very few galaxies 
with logM* > 11.3 M Q , and virtually all logM* > 11.5 M© 
galaxies are in the n > 7 bin. That is, there is a correlation 
between M* and n within our sample. Because the trend is 
towards lower values of M„/Md with increasing n, this cor- 
relation contributes to the apparent mass dependence of the 
dynamical-to-stellar mass ratio. 

Again, the dynamical mass estimates used thus far have 
been derived under the assumption of dynamical homology 
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FIG. 4. — Comparing stellar and dynamical mass estimates, accounting 
for both structural and dynamical non-homology. — The difference between 
this Figure and Figure [T] is that we have used GMd „ = Ky(ti)akT^, with 
Ky(n) defined as in Equation [5] otherwise all symbols and their meanings 
are as in Figure[T] The key point to be made from this Figure, in comparison 
to Figure [l] is that allowing for non-homology makes a big difference to the 
inferred dynamical masses, and so to the correspondence between stellar and 
dynamical masses. Further, we note that the apparent inconsistency whereby 
M, > M A for M, < 10 10 5 Mq galaxies seen in Figure [T] is no longer 
apparent. After accounting for structure-dependent differences in galaxies' 
dynamics, the relation between M„ and M& „ is much more nearly linear. 
However it remains true that the difference between stellar and dynamical 
mass appears to grow with increasing mass. 

(ie. Ky is equal to 4 for all galaxies). The solid line in Figure 
03 shows the expected shape of the relation between M*/Md 
and n, derived from Equation [3] The agreement between the 
observed relation between /M& and n and the expectations 
from this (admittedly simple) dynamical model immediately 
suggests that non-homology has an imp ortant effect on dy- 
namic al m ass estimates. We no te that both Prugniel & Simien 
( 1996) and Truii llo et al J (12004!) have made a similar argument 
for the importance of non-homology based on the fundamen- 
tal plane of elliptical galaxies (see also Section[6]below). 

4. RESULTS II— COMPARING STELLAR AND 
DYNAMICAL MASS ESTIMATES ACCOUNTING 
FOR DYNAMICAL NON-HOMOLOGY 

In this Section, we investigate the potential impact of non- 
homology on the agreement between stellar and dynamical 
mass estimates. To this end, we repeat the comparisons be- 
tween stellar and dynamical mass estimates presented above, 
using the structure corrected dynamical mass estimator, Mj „, 

in place of the simple estimate We quanti fy th e relation 
between an d Mj,,, for our sample in Section l4~Tl and show 
in Section 14.21 that allowing for non-homology considerably 
improves the correspondence between stellar and dynamical 
mass estimates. 

4. 1 . The Relation Between Stellar and Dynamical Mass 

In Figure |4j we show the relation between stellar and dy- 
namical mass for the |G09| sample, using the structure cor- 
rected dynamical mass estimator, Md.„; this Figure should be 



compared to Figure Q] It is immediately obvious that the cor- 
relation between M* and M± n is much closer to linear than 
that between M* and the simple dynamical mass, M&. Fur- 
ther, we note that the results are now logically consistent, in 
that < Md „ for almost all galaxies. This is our most ba- 
sic result: structure-dependent differences in galaxy dynamics 
can have a big impact on the inferred dynamical mass, and so 
the stellar-to-dynamical mass ratio. 

The best-fit parameters for the M^-Mj „ relation are a = 
0.92 ±0.007 and b n = -0.23 ± 0.004, where we have used 
bootstrap resampling to estimate the statistical uncertainty. 
While the statistical errors in the fit parameters are impres- 
sively small, systematic errors are sure to dominate. To see 
this, consider what would happen if we were to impose a 
stellar mass limit > My lm in Figure HI we would only in- 
clude those galaxies with Md,„ < Mn m that have high values 
of M*/Md,„; similarly, we would exclude those galaxies with 
Md.n < Mn m that have low values of M* /M±„ . This would lead 
to a significantly shallower best-fit slope to the M*/M±„ rela- 
tion. As a specific example, if we were only to consider galax- 
ies from the lG09l sample with M* > 10 10 8 M , we would find 
oc M^ f t 6 . In this context, it is worth noting both that the 
scatter in the M^-Mj,,, relation appears appears to increase 
with decreasing mass, and that the lowest mass points in Fig- 
ure|4]lie systematically above the best fit line. 

Note that, while we have phrased this in terms of a mass 
selection effect, at least part of this effect is related to how 
we have fit for the slope of the M*-M±„ relation; viz., by 
minimizing the mean perpendicular distance between the data 
and the best-fit line. If we were instead to fit by minimizing 
the mean offset in M^/M^ n as a function of M*, we would 
reduce our vulnerability to this effect. 

There are thus two effects that have the potential to sys- 
tematically bias the measured slope of the M*-M± n relation. 
We can obtain a simple estimate for the systematic error on 
the parameters a and bn by re-fitting the M*-Md,„ relations 
in different ways. If we fit by minimizing the mean vertical 
offset, AM*, from the best- fit M#/M& in relation, we find 
a = 0.85 and b\\ = -0.20. If instead we fit by minimizing 
the mean horizontal offset, AMj,,,, we find a = 1.00 and 
b\\ = -0.26. That is, the systematic uncertainties related 
to the technique used to fit the M^-Mj,,, relation (for this 
sample) are on the order Aa = 0.08 and Abu = 0.03. 

What about the systematic biases due to the particular mass 
distribution of galaxies in the G09 sample? To explore the 
importance of these effects in our measurement of the slope 
of the M*-Md,„ relation, we have tried re-fitting the M„-Md,„ 
relation, weighting each point according to its stellar mass. 
The specific weights have been derived through a compari- 
son between the mass distribution of galaxies within the IG09I 
sampl e, in bins of AM, =0.1 dex, and the z ~ mass func- 
tion of Bell et al] (120031) . We have chosen the weight for each 
gal axies so that that the weighted stellar mass distribution of 
the lG09l sample matches the 'real' stellar mass function. This 
weighting scheme is akin to 1 /V max weighting, inasmuch as if 
one were able to derive V max values for the sample, one would 
hope to obtain similar values. 

Re-fitting the G09 sample using these weights, we find 
a = 1.00 ±0.05. The larger random error on this value in 
comparison to our fiducial values stems from the fact that 
the lower mass galaxies are given much greater weights (by 
several orders of magnitude); the inclusion/exclusion of these 
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Fig. 5. — Comparing stellar and dynamical mass estimates accounting for structure-dependent differences in galaxy dynamics. — The difference between this 
Figure and Figuref2]is that we have used GMi„ = AV(«)o"o^e' with Ky(n) defined as in EquationfJ] otherwise all symbols and their meanings as is in FigurefT] 
After accounting for structure dependent differences in galaxy dynamics, the apparent trends in M^/Mj ,, with stellar mass and Sersic index are substantially 
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disappear (see also Figures [7jand|9). 
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points in the bootstrap resampling thus has a major impact on 
the best-fit slope. The fact that using these weightings gives 
an almost perfectly linear relation between M* and Mj,,, is 
striking, but it is important to remember that this fit is based 
almost entirely on the relatively small number of M* < 10 11 
Mq galaxies in the sample. In Appendix iBl we perform the 
same analysis for a more general galaxy sample, drawn from 
the NYU VAGC, which provides a much better sampling 
of the true galaxy mass function. Using the weighting 
scheme described above, the relative weights of galaxies with 
10.2 < logM*/M < 11.7 differ by a factor of only 10. For 
this galaxy sample, we find a = 0.930 ± 0.004 without 
weighting, compared to a = 0.933 ± 0.007 with weighting. 
While the data are consistent with a linear relation between 
and Md,„, they thus prefer a slightly shallower relation. For 



thelG09l sample, we find a = 0.93 ± 0.007 (±0.07), and b u = 
-0.23 ±0.004 (±0.03). This should be compared to the values 
of a = 0.73 ± 0.006 and b u = -0.14 ± 0.003 Cftos) for 

the simple dynamical mass estimate, M^. While the mass dis- 
tribution of galaxies within the IG09I sample can in principle 
induce a large bias in the measured slope of the M^-Md,,, re- 
lation, our best fit value is in fact consistent with that derived 
from a more general galaxy sample, in which these effects 
play a far smaller role. We will explore the potential role of 
other sample selection effects in Section l5T2l 

In Figure ID we show the trends in /M± n with mass, 
structure, dynamics, and density; this Figure should be com- 
pared to Figure |2] For each of these four parameters, the 
trends in M^/Mj,,, are significantly weaker than what we have 
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seen for M*/M d . The net differential trend across the sample 
is now on the order of 0.2 dex or so, as compared to 0.5 dex for 
M*/M&. While there is still a strongly statistically significant 
trend in M» /M±„ with a, the trend with surface density is now 
only significant at the 4a level. While we do still see signs of 
a trend in M*/M^„ with n, this trend is not statistically signif- 
icant, at least for the sample as a whole. We discuss this point 
further in the next Section. 

4.2. Does M*jM& n Depend on Mass, or Structure, or Both? 

In Figure |6] we return to the issue of the n- and ./In- 
dependence of M*/M d ,„; this Figure should be compared to 
Figure [3] In Figure^, we show the median relation between 
M*/M<i t „ and M* in bins of n. Again, the trends in M»/Md.„ 
with M» for the different n-bins are parallel, but offset from 
one another. In each of the 3 < n < 7 bins, we find that 
M»/Md. n scales approximately as Mj 1 ; that is, roughly con- 
sistent with the scaling that we see for the sample as a whole. 
Figure |6j? shows the median relation between M»/Md.„ as a 
function of n for different bins in M*. While the trend in 
M„/Md.„ with n is substantially weaker than we saw using 
the simple dynamical mass, Md, we still see that M*/M<i t „ 
varies with n; if anything, it would appear that by using the 
prescription for Ky(ri) given in Equation [3] we have overcor- 
rected for the effects of non-homology. Without detailed dy- 
namical modeling, however, we have no means of refining the 
model used to derive Equation [3] (We will discuss this point 
further in Section[6]) 

In other words, we have shown that accounting for struc- 
tural and dynamical homology significantly improves the 
agreement between stellar and dynamical mass estimates as 
a function of Sersic index, n, but we have not unambiguously 
shown whether or not M*/Md.„ depends on galaxy structure 
— nor can we. 

5. RESULTS III — EXPLORING POTENTIAL BIASES IN 

M»/M d ,„ 

In this Section we discuss three general classes of biases 
that may affect the results we have presented in Section |4] 
first, systematic biases in th e Sersic fits that we use to derive 
M* and M^ n (Section 15.11 see also Appendix IB!: then, the 
possibili ty o f severe selection effects for the IG09I sample 
(Section 15.21 see also Appendix IB!; and finally, systematic 
effects associat ed w ith the estimation of stellar mass-to-light 
ratios (Section 15.3b . We will show very good consistency 
between the values of M* and Md.„; this is not the case for 
the simple estimate M d . As in the previous Section, this 
consistency provides strong circumstantial evidence — but 
not proof beyond a reasonable doubt — that there are no 
significant biases in either measurement. 

5.1. Looking for Possible Observational Biases 

We explore the possibility of serious observational biases 
in Figure |7] In each panel of this Figure, we plot M»/Md.„ 
as a function of a basic observable: namely (left to right), 
apparent magnitude, apparent size, and redshift. The lines 
in each panel show the median relation for the same bins of 
stellar mass as are shown in Figures |3j) and|6j). 

By a similar argument to the one given in Section l4~2l if the 
apparent mass dependence of the ratio M*/M^ n were driven 
by observational effects that are strong functions of apparent 
brightness or size, then we would expect there to be clear 



trends in M*/M d ,„ for each of the stellar mass bins in Figure 
[7] This is not obviously the case. Fitting to the data in each 
stellar mass bin, the gradients of the best-fit relation between 
M„/Md.„ and all of apparent magnitude, apparent size, and 
redshift are consistent with zero; this is true for each stellar 
mass bin individually (typically within la), as well as for the 
sample as a whole (within 1.2-1. 4a). 

We note that the same is not true using, for example, the 
Sersic structural parameters given in the NYU VAGC to 
derive M» and M dj „ (as we do in Appendix IB1. In this case, 
we do see a weakly statistically significant gradient (at the 
level of 4a) with observed size. That is, this kind of test is 
indeed able to (weakly) detect mild systematic errors in the 
Sersic-fit parameters on the order of 10-20 %. The results in 
Figure|7]thus argue against the idea that there are any serious 
biases affecting the measurement of M„ or M d ,„ (or, more 
accurately, the ratio M»/Md.„) associated with the Sersic-fit 
structural parameters used to derive these values. 

As can be seen in Figure [5] there is a statistically signifi- 
cant correlation between M*/M di „ and velocity dispersion, a. 
Using the same argument as above, it is conceivable that this 
could be produced by a systematic bias in the measured val- 
ues of a. The observed trend of ~ -0.12 dex inM*/Md ; „ over 
~ 0.5 dex in a could be entirely explained by a ~ 0.06 dex 
differential bias in the measured values of a. In connection 
with this point, we stress that a com parison between th e SDSS 
measured values of a and those of Faber et al. ( 1989) shows 
no systematic biases. This should give some confidence that 
the trend in M*/Md.„ with a is real, and not a product of ob- 
servational biases. 

5.2. Looking for Potential Sample Selection Effects 

As we have described in Section 12.41 the IG09I sample has 
been constructed in such a way that massive galaxies in gen- 
eral, and in particular central galaxies in very massive ha- 
los, are drastically overrepresented in comparison to the gen- 
eral field population. If there are systematic differences in 
M*/M^„ as a function of, for example, environment or star 
formation activity, there is thus a very real danger that sample 
selection effects may play an important role in shaping our 
results. 

We expl ore t his issue in Figure |8] in which we have 
divided the G09 sample in central/satellite and non-/emission 
subsamples, in order to look for differences between these 
populations. Here, we have selected 'non-emission' galaxies 
as Ha and H/3 in absorption. The 'emission' subsample 
contains both star forming galaxies and AGN hosts. 

In the upper panels of Figure [8] it is clear that we find 
essentially identical relations between M* and M d ,„ for the 
central and satellite galaxy subsamples. Given that, as we 
have shown in Figure [6] M*/Md.„ varies with both M» and n, 
it is not all that surprising that central and satellite galaxies 
show the same M*-M d ,„ relation: not only have the two 
subsamples have been constructed to be matched in . IG09I 
have shown that there are no structural differences between 
satellites and centrals at fixed mass. 

However, we do find very slightly different M*-Md,„ 
relations for the emission and non-emission subsamples. For- 
mally, the two relations are slightly offset from one another, 
at the level of 0.05 dex; the gradients of the two relations 
differ at the 2a level. That said, for the 10.5 < logM* < 12 
range spanned by the sample, the difference between the 
best-fit relations for each subsample and that for the sample 
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FIG. 7. — Signs of observational biases? — Each panel shows M* /M±„ a function of a direct observable. Within each panel, the solid lines show the median 
relation in bins of stellar mass; these bins are the same as those shown in the right-hand panels of Figuresf3]and|6] The points show the data themselves; points 
that fall outside the range of each panel are shown as a small grey plus. Within the G09 sample, there are correlations between Sersic index and each of the 
observed quantities shown in this Figure; however, there are only very weak trends in M t /Mj „ with any of these observables. Moreover, each of the stellar 
mass bins follows essentially the same median relation. This suggests that neither the stellar nor dynamical mass estimates are obviously seriously biased by 
systematic errors in the Sersic fits. 



as a whole are < 0.03 dex. We also note that the apparent 
offset disappears (at least within 2a) if we consider only the 
n > 2 galaxies in both subsamples. Moreover, we do not 
find any differences in M*/Md.„ at fixed M* and n between 
these different subsamples — within statistical uncertainties, 
each subsample shows the same behavior as is seen in Figure 
[6] This suggests that the apparent offset between for the 
emission and non-emission galaxies is driven by the different 
distributions of Sersic indices within the emission and non- 
emission subsamples, rather than any difference in M JMa „ 
intrinsically related to galaxy activity. (See also Section l531 l. 

That is, our results do not appear to be strongly influenced 
by the relative numbers of central/satellite galaxies or of non- 
/emission galaxies in our sample. The above analysis does 
suggest, however, that the measured relation between M* 
and M& „ is sensitive to the joint M„-n distribution within the 
sample. This is a direct consequence of the fact that M*/M<i t „ 
depends on both mass and Sersic index (Figure [§}. We have 
conside red biases associated with the mass distribution in 
Section 14.11 But note that if, at fixed mass, the distribution 
of n within the IG09I sample differs significantly from the 
'true' distribution for a general field population, then we may 
therefore find a very different slope for the M*/Md,„ relation. 

For this reason, in Appendix |B] we repeat our M»-Md „ 
comparisons for a more general sample of 0.035 < z < 0.08 
field galaxies, using the structural parameters given in the 
NYU VAGC0 The best fit logarithmic slope of the M*-M d „ 
relation for 0.035 < z < 0.08 field gal axies i s a = 0.91 ±0.003, 
compared to a = .92 ± 0.009 for the lG09l sample. This sug- 
gests that the IG09I sample is not grossly biased in terms of the 

8 As we have already remarked, the VAGC Sersic fits have been shown 
to suffer from systematic errors arising from background oversubtraction. 
But, as we also show in Appendix [E] it turns out that the ratio M*/M t i „ is 
extremely robust to random or systematic errors in the fitting of structural 
parameters, so long as the apparent magnitude, effective radius, and Sersic 
index are derived consistently, and the term Ky(n) is included in the defini- 
tion of the dynamical mass (See also Section|6j. Further, we find very similar 
M*-M,i „ relations for the G09 sample using either VAGC or the GKj values 
for the structural par amet ers. Any large differences in the measured M*-M,j „ 
relation between the G09 sample and the general field sample would therefore 
necessarily be a product of selection effects. 



distribution of n at fixed mass. (Here again, it is significant 
that IG09I have found that, at fixed mass, there are no struc- 
tural differences between satellite and central galaxies.) 

To summarize the results of this section, then, separate anal- 
ysis of central/ satellite and non-/emission galaxies suggest 
that our results are not strongly affected by selection effects 
associated with these properties. Furthermore, although the 
observed slope of the M*/M<i t „ relation is in principle sensi- 
tive to the joint M*-n distribution within the sample, we find 
very little difference between the IG09I sample and a more gen- 
eral field galaxy sample. Finally, we stress that we do not find 
any evidence that selection effects have an important impact 
on the results shown in Figure |6j ie., the observation that the 
ratio M*/Md „ depends on both n (at fixed M*) and on M* (at 
fixed n). 

5.3. Looking for Biases in the Stellar Mass-to-Light Ratio 
Estimates 

As we have stated in the Introduction, our primary moti- 
vation for comparing stellar and dynamical mass estimates is 
to validate the stellar mass estimates. We explore this issue 
in Figure |9] In each panel of this Figure, we plot the ratio 
M*/M(j,« as a function of a different property of the stellar 
population. The solid lines in each panel show the median re- 
lation for the same stellar mass bins shown in Figures [3] and 
[6] The upper panels plot M»/Md.„ as a function of a direct 
observable; the lower panels plot M*/Md : „ as a function of 
a derived property. Note that the stellar mass estimates we 
have used were derived from the ugriz photometry, rather than 
spectra. The measured values of M*/Md.„ are thus formally 
independent of the three spectral measurements shown in the 
top panels. F urther, note that the age estimate that we show 
is taken from Kauff mann et al.l d2003al) ; these values are also 
derived from the spectra. The dust obscuration and specific 
star formation rates (SSFRs) shown come from the SED fits 
used to derive M»/L; these values are thus self-consistently 
derived. 

We have color-coded the data in Figure [9]according to thei r 
spectral classification as given by Brinchma nn et al.l (2004); 
viz.: star forming (blue); low S:N star forming (cyan); com- 
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FIG. 8. — Possible sample selection effects? — The G09 sample is not necessarily representative of the general galaxy population, in terms of the distribution 
of masses, environments, or_star formation activity. To explore the potential role of these effects, each panel of this Figure shows the M*-M,i „ comparison 
for different subsets of the G09 sample, distinguishing between central and satellite galaxies (upper panels) and between spectroscopically emission and non- 
emissi on galaxies (lower panels). For the lower panels, the spectroscopic classification is based on the BPT diagram, following the scheme of Brinchmann et al. 
12004); t he 'emission' sample inc ludes both AGN-dominated and composite spectra galaxies. In the upper panels, we do not show those satellite/central that do 
not have Brinchmann et al. 12004) spectral classifications; that is, the same samples are plotted in both the upper and the lower panels. The dashed grey line in 
each panel shows the best fit relation for the entire G09. We find a consistent M*—M&j, relation for all four of these subsamples, as well as the G09 as a whole. 
While there is possibly a slight offset in M*/Ma,„ between star-forming and passive galaxies, we note that this disappears if we only consider n > 2 galaxies. 
That is, this offset appears to be due to the different n-distributions of the star-forming and passive samples, rather than an intrinsic difference between the values 
of M t /Mi j, for emission and non-emission galaxies (see also Figure [§}. We also note that for each of the subsamples shown, within statistical errors, we find 
consistent behavior in M t jM& „ at fixed M„ and n as is shown in Figure|6] We therefore conclude that selection effects do not play a major role in shaping our 
results (see also Appendix IB"V 



posite (yellow), AGN (orange), and non-emission (red). In 
general, we see little if any differences in the values of 
M»/Md „ for different stellar populations. There are not obvi- 
ously large differences between the mean values of M»/Md.„ 
for AGN hosts, star forming galaxies, or non-emission galax- 
ies. 

Within each panel we give the gradient of the best-fit line 
for the whole sample. These values are all statistically con- 
sistent with zero. The only possible exceptions to this rule are 
the gradients in M» /M^ n as a function of Ha equivalent width 



(EW) and as a function of M*/L r , both of which are non-zero 
at the ~ 3cr level. Just as there is little if any trend for the 
sample as a whole, there are no statistically significant trends 
for any of the individual mass bins. 

We can quantify the degree of correspondence between stel- 
lar and dynamical mass estimates by considering the differen- 
tial bias between galaxies over the range of each of the prop- 
erties shown in Figure [9] Looking at the median relations 
shown for each mass bin suggests that the magnitude of such 
differential biases are at most 0.2 dex. We can obtain similar 
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FIG. 9. — Stellar population-dependent effects? — Each panel plots M*/Mi„ as a function of a stellar population diagnostic. The upper panels show directly 
observed quantities; the lower panels show stellar population parameters derived from synthetic stellar population modeling. In these panels, the mean stellar 
age shown is that given by Kauffmann et al. 1 2003a), which are based on spectra, and are thus formally indepe nden t of M*/L. The dust extinction, specific 
star formation rates, and M„/Ls are all self-consistently derived from the same SED fits, as described in Section [231 Within each panel, individual points are 
color-coded according to spectral classifications; specifically: star forming (blue), low S/N star forming (cyan), composite (yellow), AGN-dominated (orange), 
and 'unclassifiable' (ie. non-emission; red) galaxies. Within each panel, the black points with error bars show the median and 16/84 percentiles of M^/Mj „ in 
bins. We see no statistically significant systematic differences in M* /Mi „ for galaxies with different stellar populations or star formation histories. This argues 
against there being major problems with the stellar population models used to estimate M*/L. At 99 % confidence, these results suggest that any differential 
biases as a function of the parameters shown is < 0.12 dex (rj 40 %). 



estimates for the sample as a whole using the fit parameters 
given in each panel of Figure|9] Taking the 3<r statistical lim- 
its on the slopes of these relations, we find that the differential 
effects across the full range of the sample are < 0.1 dex for 
D„(4000), RS EW, and restframe color; and < 0.15 dex for 
Ha EW, age, dust extinction, SSFR, and M*/L. 

The same is not true using the simple estimate of dynam- 
ical mass, Ma, in place of Md,„. As might be expected from 

comparing Figures Q] and [4] we find that the M*/Md relations 
for different mass bins are largely parallel, but significantly 
offset from one another. Then, because of correlations 
between mass and activity, we also find significant gradients 
in the M*/Md relations for the sample as a whole; typically 
at the 10-20cr level. The size of differential biases as a 
function of all of D„(4000), RS EW, Ha EW, age, and SSFR 
are on the order 0.2-0.4 dex. We note in particular that the 
relatively strong gradie n t in M a /M* with Ha EW noticed by 
Drory, Bender & Hopp (2004) disappears when we use the 
structure corrected dynamical mass estimator, Md.„, in place 
of the simple estimate this apparent bias seems to be 
more closely linked to structure than to Ha emission per se. 
We thus find a very good correspondence between M* and 
Ma,n for galaxies in the lG09l sample, but only provided we 
account for structural and dynamical non-homology. 

5.4. Color-M J L Relations 

iBell & de Jongl (1200 lb have shown that M*/L and color 
are strongly correlated in both the optical and the NIR. This 
implies that a single color is enough to make a reason- 



able estimate of M*/L (see also, e.g., lLong hetti & S araccol 
I2009t iGallazzi & BdH 120091) . These kinds of color relations 
have since been w idely used for high redshift studies. In 
iTaylor etafl (120091) . we have shown that the values of M^/L 
that we use here correlate very strongly with A (g-i) color; 
viz.: 

logM„/L, = -0.82 + 0.83 x 0A (g-i) . (4) 

(In this expression, L,- should be understood as referring to 
the absolute luminosity in the I /-band filter; that is, the 
/-band filter redshifted to z = 0.1. The absolute magnitude 
of the sun in the '/-band is 4.58). The scatter around this 
relation is just 0.03 dex. If w e use this relation to predict 
M r /L for galaxies in the IG09I sample, we again find very 
good correspondence between M„ and Md,„ — in fact, the 
rms scatter in M»/Md.„ is unchanged. 

That is, at least from the point of view of consistency be- 
tween M„ and M^ n , it would seem that M»/Ls estimated on 
the basis of a single color are not significantly worse than es- 
timates based on full SED fits. This is significant because the 
SEDs that were used to derive the M»/Ls were corrected for 
emission lines using the SDSS spectra. This kind of correc- 
tion is not practicable for, for example, high redshift studies. 
Further, we have repeated our analysis using th e M*/L s de- 
rived from the SDSS spectra by Kauffmann et al. (2003a), and 
find similarly good agreement between M* and Md,„: at least 
on average, it would appear that color-derived M*/Ls are just 
as good as those deri ved from optical spectroscopy (see also 
Gall azzi & Bell2009h . 
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As can be seen in Figure [9] we find no statistically signif- 
icant trend in M*/Md,„ with 0A (g — i) color. Given the exis- 
tence of a relation between M*/L and color, this fact implies 
a relation between M^„/L and color. Using the color-M^/L, 
relation given above, the results shown in Figure[9]imply that: 

logM d .„/Li = (-0.54 ± 0.07) 

+ (0.80 ± 0.05) x 01 (g-i) . (5) 

(Again, L, in this expression should be understood as referring 
to the 01 i-band.) For the galaxies in our sample, the scatter 
around this relation is 0. 14 dex, almost exactly the same as the 
scatter around the M*-M±„ relation. That is, using only the 
01 (g — i) color, it is possible to predict the dynamical masses 
of the galaxies in our sample with a relative uncertainty of 
« 40%. 

6. DISCUSSION 

6.1. (Non)-Homology and Dynamical Mass Estimation 

In retrospect, it is perhaps unsurprising that is not a good 
estimate of dynamical mass. By using a Sersic-fit size and 
total magnitude measurements we have allowed for structural 
non-homology in our definition of Mj, while at the same time, 
by using k = 4 for all galaxies, we have assumed dynamical 
homology. In this sense, the definition of is not internally 
consistent. 

On the other hand, it is remarkable that the SDSS model 
sizes can be used to obtain a reasonably good dynamical mass 
estimate, under the assumption of both structural and dynami- 
cal homology. As we argue in AppendixlBl the n-dependence 
of Ky(n) and cthe ovariance between n and the combination 
R e /L are very nearly equal and opposite to one another. 
This leaves M± n /L oc Ky(n)a 2 R e /L remarkably insensitive to 
errors in the size and total magnitude measurements, provided 
that 1.) m tot , O e , and n are self-consistently derived, and 2.) 
Md „ is derived using the appropriate value of n to calculate 
the dynamics term k or Ky(ri). This means that, for example, 
assuming a De Vaucouleurs profile (ie., n = 4) to derive m T and 
9 e for all galaxies ensures that M&/L w Md,„/L, provided that 
the appropriate value of Ky is used; viz., k « Ky(n = 4) = 4.62. 
This point is amply demonstrated in Figure QT| That is, 
M^ n /L appears to be insensitive to errors or assumptions in n, 
so long as everything is consistently derived self-consistently. 

6.2. Ky(ri) and Dynamical Mass Estimation 

The inclusion of the term Ky(n) makes the structure- 
corrected dynamical mass estimator M^„ explicitly model de- 
pendentQ Further, this prescription for Ky(n) has been de- 
rived under very simple and idealized assumptions {viz. a sin- 
gle component, spherical, and dynamically isotropic distribu- 
tion), and so can only be regarded as approximate. 

That said, more sophisticated dynamical models can give 
an indication as to how large these effects might be. For 
the case of anisotropy, the effects on the value of Ky(n) 
are on the order of < 0. 1 dex, and become le ss impor- 
tant for larger n (see, e.g., ICiotti & Lan zoni 2001). Further, 
iBertin. Ciotti & Del Principd (120021) argue that the galaxy dy- 
namics (or, more accurately, the value of Ky) close to the 

9 By the same token, the simple dynamical mass estimator is also 
model dependent, inasmuch as it assumes homology, which is patently 
wrong. 



galaxy center are in principle rather sensitive to the precise 
shape of the total mass distribution. Their results suggest that 
this effect is on the order of < 0.1 dex (see their Figure D.l). 

In order to probe the dependence of our conclusions on 
the assumed form of Ky(ri), we have also trial ed using an 
alterna te prescription for Ky(n), given by Cappellari et al. 
(2006, see their Equation 20). The main difference between 
this prescription and the one given in Equation [3] is that it is 
phrased in terms of the observed velocity dispersion within 
the effective radius, cr e , rather than the central velocity dis- 
persion (To. This prescription thus has a different dependence 
on dynamical isotropy and the dark matter profile, and so 
provides an indirect m eans of probing the i mportance of 
these effects. Using the Cap pellari et all (12006b prescription, 
we find qualitatively and quantitatively similar results: we 
find that M* oc M^ 5 ; we still see that M»/M d .„ depends on 
M» at fixed n, and on n at fixed M„ (although this dependence 
is somewhat shallower); we see no trends in M*/Md.„ 
with apparent magnitude or redshift, but a weak trend with 
observed size; and, at fixed mass, we see no statistically 
significant trends in M„/L with stellar population parameters. 

That is, while it is virtually certain that the model used to 
derive the prescription for Ky(n) given in Equation|3]is wrong 
in several important respects, it seems unlikely that account- 
ing for those effects that are ignored in the model would have 
a drastic effect on our results and conclusions. As we have re- 
peatedly stressed, detailed dynamical modeling is necessary 
better constrain the 'true' values of Ky for individual galax- 
ies. 

6.3. Comparison to Previous Studies 
6.3.1. Comparison to Other SDSS Studies 

In comparison to other studies of the relation between 
stellar and dynamical mass estimates based on SDSS data, 
we find con siderably less variation i n M*/ M±„ with mass. In 
the case of iDrorv. Bender & Hoppl (120041) . this difference is 
simply due to the fact that we account for non-homology in 
the derivation of Md,„; using the simple mass estimate Mj, 
we have verifie d that we are able to reproduce their results. 
In the case of Gall azzi et alJ (l2006h . there is the additional 
complication that they use qualitatively different measures of 
total flux and size: specifically, the Petrosian magnitudes and 
half-light radii given in the basic SDSS catalog, which are 
derived directly from the observed curves of growth. Again, 
we have verified that we can reproduce their results using the 
s ame measurements. 

iGallazzi et alJ ( 12006b also split their galaxy sample into 
bins of Sersic index, and find similar slopes to the M*-Md 
relation for each subsample; they find that the logarithmic 
slope of the relation varies from 0.847 to 0.801 between n = 3 
and n = 5.5. The size of this variation is entirely consistent 
with the results we have shown in Figure [3] They use this 
fact to argue that non-homology does not have a significant 
impact on the slope of the global M*-Md,„ relation. However, 
as we have also shown in Figure[3] while the M^-M^ relations 
for each bin in n are parallel, they are significantly offset 
from one another. It is this through this offset, combined with 
a correlation between and n, that non-homology affects 
the slope of the M*-M± n relation; Gall azzi et alJ (|2006) make 
no mention o f such an offset. Wh ile our conclusions differ 
with those of Galla zzi et alJ (120061) . our results are thus not 
obviously inconsistent. 
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6.3.2. Comparison to Detailed Dynamical Modeling Results 

Cappella ri et al.1 (|2006) have also argued against the idea 
that non-homology has an important impact on dynamical 
mass estimates. This argument was based on dynamical mass- 
to-light ratios derived from detailed 2D and 3D modeling 
25 str ucturally early t ype g alaxies from the SAURON sam- 
ple. ICappellari et al.l (120061) compared the dynamical mass- 
to-light ratios, so derived, to the simple virial mass estima- 
tor Md/L oc CT 2 /?DeV /Loev, where Ru e v and Lo e v were derived 
from De Vaucouleurs profile fits. They found no evidence for 
an n-dependent offset between these tw o quantities. 

The fact that Cappellari et al.l d2006 i) use De Vaucouleurs- 
fit sizes and magnitudes is significant: as we have argued 
above and in Appendix IB~1 when using De Vaucouleurs fits 
to derive size and magnitude measurements, it is appropri- 
ate to use a constant k in the definition of Mj. That is, by 
using De Vaucouleurs- fit sizes and masses to define (M^/L), 
Capp ellari et al.l (120061) effectively guaranteed, almost by con- 
struction, that they would find no structure-dependence when 
comparing the 'simple' dynamical mass-to-light ratio to that 
derived from full dynamical modeling. 

6.3.3. Comparison to Fundamental Plane Studies 

As we have mentioned in the Introduction, the fundamental 
plane can be thought of as measuring the variation in the dy- 
namical mass-to-light ratios of early type galaxies as a func- 
tion of velocity dispersion, luminosity, or mass (see, e.g . , 
iDressler et ail 119871: lJ0rgensen, Franx & Kiaergaardl fl996[) . 
How do our derived values of Md.„ compare to those de- 
rived from the fundamental plane? To address this que stion, 
we selected the non-emission galaxies from within the IG09I 
sample with n > 2.5. For these galaxies, we find M^ n jL oc 
0_o.88io.O6. For comparison, Jtorge nsen. Fran x & Kiaergaard 
(1996) find o- 2 R DeV /L 

Dev oc cr , where again, /?DeV and Lo e v 
have been derived via De Vaucouleurs fits. Again, the co- 
variance between the fit values of m r , 9 e , and n mean that 
c 2 /?DeV /Loev ~ M^„/L. Our structure-corrected dynamical 
mass estimates are thus in good agreement with those derived 
fro m the fundamental plan e. (Parenthetically, we also note 
that ICappellari e t al. (2006) found that their dynamical mass- 
to-light ratios derived from detailed dynamical fits scaled as 
a 82 , which is also consistent with our results.) 

We have also considered how the dynamical-to-stellar mass 
ratio, Md^/M* varies with a for this same sample of early 
type galaxies: we find M d ,„/M» oc a °- 50 ± 006 This would 
suggest that less than half of the tilt of the fundamental plane 
is due to variations in the mass-to-light ratios of early type 
galax i es as a function of a (cf.. e.g.. IPrugniel & S imien 
Il99l iTrurliio et al.ll2004t lAllanson et al.ll2009l) . We present 
this result only for completeness; proper interpretation of this 
result requires much more detailed analysis, and is beyond 
the scope of this work. 

We also note th at both IPrugniel & Simien! d 19961) and 
Truiill oet al.l (|2004) have made a very similar argument for 
the importance of non-homology in estimating dynamical 
masses as we have made in Sections[3]and|4] based on dynam- 
ical mass-to-light ratios derived from the fundamental plane. 
Our analysis based on the correspondence between M* and 
Md „ is complementary to theirs in two ways. First, their anal- 
yses were specific to early type galaxies; we have thus ex- 
tended their result to the general galaxy population. Secondly, 
both authors focussed on M^/L, rather than M»/Md.„; that 



is, neither of these authors considered the relation between 
galaxies' stellar and dynamical masses. 

6.4. Interpretation 

Turning now to the interpretation of our results, the remark- 
able consistency between stellar and dynamical mass esti- 
mates shows two things. First, it strongly suggests that the 
measurements of M* and Md,„ are both meaningful and rela- 
tively robust. In particular, our results indicate that it is possi- 
ble to derive stellar mass estimates without strong differential 
biases as a function of age, dust, SSFR, or M*/L, based only 
on broadband optical photometry (or indeed on a single opti- 
cal color). 

Secondly, it implies that intrinsic variations in the stellar- 
to-dynamical mass ratio (due to, e.g., variations in the dark- 
to-stellar mass ratio, variations in the IMF, or dynamical dif- 
ferences beyond the simple non-homology considered here) 
as a function of stellar mass, galaxy structure, and star for- 
mation rate/history are either small, or conspire to leave the 
inferred values of M^jM^„ relatively unchanged. The rest of 
this section will be devoted to discussion of this result. 

Since Md.„ is an estimate of total mass, it can only be in- 
terpreted as an upper bound on the true stellar mass. Because 
we have no a priori means of separating out the relative con- 
tributions of luminous and non-luminous mass to M^ n , the 
relation between M» and Md.„ is complicated by degeneracies 
between the relative contributions of gas and dark matter, as 
well as uncertainties in the low mass shape of the IMF. 

The simplest way to interpret the non-linearity of the 
M»-Md.„ relation is as indicating a greater central dark matter 
fraction for higher mass galaxies, in qualitative agreement 
with theoretical expectations. Using simple argu ments based 
on t he observed dyna mics of elliptical galaxies. iFranxl (119931) 
and lKochanekl £1994) have argued that accounting for a dark 
matter halo implies |logM*/M d .„| ~ 0.14-0.18 dex. This 
would go a long way towards explaining the ~ 0.23 ± 0.03 
offset that we have observed. 

We can es t imate gas masses using the prescription given by 
Zhang et al.l (120091) . These authors have used a sample of rel- 
atively low mass SDSS galaxies with literature Hi masses to 
derive a prescription for Mm / M» as a function of (g - r) color 
and stellar surface density. Using this prescription to derive 
baryonic mass estimates, Mt, ar = M» +Mhi, reduces the size of 
the offset between Mb al and Md,„ by 0.05 dex to -0. 1 8 dex, and 
brings the logari thmic slope of the M \, m /M± n relation to 0.95. 
The fact that the Zhang et al. (2009) relation has been derived 
for very different galaxies to the ones we consider here means 
that this result should be interpreted with caution. Even so, 
it is striking that, taken together, the estimated contributions 
of Hi and dark matter almost perfectly explain the observed 
offset between M* and Md,„, and imply only a mild trend in 
M bal -/M d ,„ with mass: M bal -/M d ,„ oc M^f 5 . 

Then there is th e matter of the IMF. The e ffect of adopting 
a ISalpeteri (119551) IMF rather than that of IChabrieri (120031) 
would be approximately to scale all our values of M* up by 
0.22 dex. For a linear M*-Md,„ relation (which our data are 
only marginally consistent with), this would leave virtually 
no room for dark matter or gas in the centers of galaxies in 
our sample. For the slightly less-than-linear relation preferred 
by our data, this would imply that > M^„ for galaxies 
withM, < 10 11 M Q , which is logically inconsistent. Thus we 
can sa y that, at best, o ur results are only marginally consistent 
with a Salpeter (1955) IMF. Accounting for dark matter, our 
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results are also weakly inconsiste nt with a ' diet S alp eter' IMF, 
and co mpletely consistent with a Kroupa (2003) or [Chabrier 
(2003) IMF. (Parenthetically, we also note that the results 
shown in Figure|9]can also provide a weak constraint on vari- 
ations in the IMF as a function of star formation rate/history.) 

Finally, we note that the observed scatter around the 
M„-Md.„ relation is rather small: just 0.13 dex. We argue in 
Appendix IB1 that the ratio M*/Md.„ oc a 2 R e /L is remarkably 
insensitive to errors in the Sersic-fit parameters, provided 
that they are consistently derived, and that dynamical 
non-homology is taken into account. This implies that the 
uncertainties in M„jM^ n are dominated by errors in the 
measurement of <jq and M»/L. The mean formal uncertainty 
in a for our galaxy sample is 0.034 dex. We estimate the 
mean random error in M* /L to be on the order of 0. 1 dex; this 
is the random scatter between the SED-fit M*/Ls used here 
and the spectrally derived M r /Ls given by Kauffman n et al.l 
(2003a)H Adding these errors in quadrature (ie. neglecting 
correlations between a and M*/L at fixed M») produces an 
uncertainty in M*/M& in of 0.12 dex. This would imply that, 
at fixed M*, the intrinsic scatter in M*/Md,„ is potentially 
very small indeed: < 0.04 dex. 

7. SUMMARY 

The central focus of this work has been the degree of con- 
sistency between stellar and dynamical mass estimates, based 
on the latest generation of data products from the SDSS. We 
have shown that structural differences in galaxy dynamics 
can have a large impact on the estimated values of dynamical 
mass, and so on the degree of correspondence between stellar 
and dynamical mass (Section |3j Figures Q] and @). Provided 
we account for structure-dependent differences in galaxy dy- 
namics (using the term Ky(n), as defined in Equations |2] and 
H), we find very good agreement between the inferred stellar 
and dynamical masses of galaxies within the G09 sample. 

Our analysis is based on the carefully-constructed satel- 
lite/central galaxy sample of G09, making selection effects 
a major potential concern. However, we find no signs of 
major differences in the relation between M» and M^ n for 
central/satellite or non-/emission galaxies within the sample, 
suggesting that our res ults are not seriously affected by 
selection effects (Section [5.21 Figure |8). Moreover, we find 
qualitatively and quantitatively similar results analyzing a 
more general sample of 0.035 < z < 0.08 galaxies, using the 
best-fit Sersic parameters given in the NYU VAGC, or using 
the De Vaucouleurs/exponential model fit parameters given 
in the basic SDSS catalog (AppendixIBl. 

We fin d tha t th e rat io M»/Md.„ varies with both M* and n 
(Sections 13.21 and |4.2| Figures [3] and [S). While the apparent 
n-dependence of M*/Md, n is sensitive to the assumed form of 
Ky{ri), changing Ky(n) cannot affect the result that M»/Md.„ 
varies with M* at fixed n. Without spatially resolved dy- 
namical information for individual galaxies, however, we 
cannot determine whether the apparent mass-dependence of 
M*/Md t „ is caused by some mass-dependent difference in 
galaxy dynamics, rather than a genuine physical difference in 
the stellar-to-dynamical mass ratios of galaxies with different 
masses. 

10 Note that this is almost certainly an underestimate of the 'true' random 
uncertainty in M* jh. Including NIR data (where the stellar population mod- 
els are the most uncertain), and properly accounting for propagation of un- 
certainties in stellar population models and the IMF, Conroy, Gunn & White 
12009) argue that the uncertainties in M*/L are on the order of 0.3 dex. 



Similarly, while we have shown very good agreement be- 
tween stellar and dynamical masses for SDSS galaxies, we 
cannot unambiguously prove that neither of these quantities 
suffers from systematic biases. On the other hand, using the 
Sersic-fit parameters given by G09, we do not see any system- 
atic variation in M*/Md.„ with observed properties like appar- 
ent magnitude, apparent size, or redshift (Section [5T1 Figure 
|7). This is not true if we use the Sersic-fit parameters given 
in the NYU VAGC, which has been shown to suffer from 
systematic errors arising due to background over-subtraction. 
That is, we have the ability to detect these sorts of errors, and 
do not see evidence for such errors for our sample. Further, 
we do not see any signs of variation in M*/M& fl for galaxies 
with different stellar populations, or for galaxies in different 
states of activity {ie., AGN hosts, star forming galaxies, or 
non-emission galaxies). 

These results, together with the good general agreement 
between M* and Md.„ provide strong circumstantial evidence 
(but not proof beyond a reasonable doubt) that there are no 
serious systematic biases in the values of M* and Md.„ that we 
use here. This implies that the assumption of non-homology 
gives the wrong dynamical mass. Further, this suggests that 
there are not strong biases in the M»/Ls we have used here: 
at 99 % confidence, the consistency between M* and M± n 
implies any differential biases in the estimate of M^/L across 
a wide range of stellar populations are at the level of < 0.12 
dex (« 40%). 
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FIG. 10. — Comp arison between the two different velocity dispersions measurements given in the SDSS DR7 catalog and those given by Faber et al. 1 1989). 
— Bernardi 1 2007) have shown that the DR5 SDSS velocity dispersions suffered from systematic biases in comparison to those from the ENEAR sample as well 
as earlier SDSS releases. For this reason, the algorithms for estimating velocity dispersions from SDSS spectra were substantially revised for DR6 and later. In 
each panel of this Figure, we compare one of the two SDSS velocity dispersion measurements to those in the 'seven samurai' catalog i Faber et al. 1989). There 
are no signs of any systematic problems with the DR7 SDSS velocity dispersions. (Note however that the scatter in these comparisons is significantly higher than 
would be expected from the formal measurement uncertainties, which are on the order of 3.5 km/s.) 



APPENDIX 

VALIDATING THE SDSS DR7 
VELOCITY DISPERSION MEASUREMENTS 

Bernardi (2007) showed that there was an inconsistency between the a-L relations for early type galaxies derived using the 
early data release (EDR) and DR5 SDSS catalogs. Further, she was able to show that the cause for this discrepancy was systematic 
biases in the DR5 velocity dispersions: in comparison to literature values from HyperLeda, the DR5 measurements slightly but 
systematic ally over-estim ated the velocity dispersions of intrinsically low-cr galaxies. For DR6+, partially in response to the 
findings of Bernardi (2007), the SDSS velocity dispersion pipelin es were substantially rev ised. The new dispersions have been 
shown to agree well with the EDR ve locity dispersio ns used by iBernardi et al.1 (l2003albl) . and thus, by implication, with the 
improved estimates for DR5 derived by Bernardi (2007)0 

In this Appendix, in order to validate the DR7 velocit y dispersions, we p resent a comparison between the velocity dispersions 
given in the basic SDSS DR7 catalog to those given by IFaber et al.l (1 1989b for elliptical galaxies in their sample. The results of 
this comparison are shown in Figure[l0] The left panel of this Figure shows the comparison for the Princeton or SpecBS values 
of <t; the right panel shows that for the Chicago or spectrold values of a. Note that the Chicago algorithm only outputs values of 
a for those galaxies that are spectroscopically classified as being early type; the three IFaber et al.l (119891) galaxies at the bottom 
of the right-hand panel are not classified as being early type, and so are not given Chicag o velocity dispersio ns. 

Within both panels, we give the median and rms difference between the SDSS and Fab er et all ([T989) velocity dispersion 
measurements. It is clear from t his Fig ure that neither of the DR7 velocity dispersions suffers from serious systematic biases 
in comparison to the Fab er et all ([1989) measurements. We note, however, that the rms scatter, which is on the order of 19 km/s, 
is considerably higher than the median formal measurement uncertainty given in the SDSS catalog, which is on the order of 
3.5 km/s for the galaxies shown in Figure [10] That is, it seems that the formal uncertainties on the SDSS velocity dispersions 
significantly underestimates the true error, at least for the relatively bright galaxies shown here. 



See http://www.sdss.org/dr7/algorithms/veldisp.html 
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FIG. 11. — Comparing dynamical and stellar mass for a general galaxy sample. — For this Figure, we have selected 0.035 < z < 0.08 galaxies with 
sciencePr imary spectra, cr b > 75 km/s, and AtT b/(T (, < 0.1. For each panel, we have derived Mi„ using either Sersic structural parameters from the 
NYU VAGC (Blanton et al. 2005a, left panel), or using De Vaucouleurs/exponential model structural parameters from the basic SDSS catalog (right panel). As 
in other Figures, the soli d line s show fits to the data; the points with error bars show the median relations in bins. The grey dashed line shows the M*-M<| „ 
relation we derive for the G09 catalog, using their Sersic structural parameters. The general M»-M,j „ relation for field galaxies is very similar to the one we find 
for the heavily-selected G09 sample. 



SELECTION EFFECTS AND SYSTEMATIC BIASES: REPEATING OUR ANALYSIS FOR A GENERAL 

GALAXY SAMPLE 

As we have repeatedly stressed in the main text, the IG09I samp le that we analyze in the main text is heavily selected. In 
order to make sure that our conclusions are not unique to the |G09 sample, in this Appendix we repeat our analysis for a more 
general galaxy sample. For this exercise, we have selected mp eU . < 17.5 galaxies with sciencePrimary spectra in the range 
0.035 < z < 0.08. As in our main analysis, we also require that the relative error on the velocity dispersion is less than 10 %, 
and that er b > 75 km/s; these selections effectively limit the sample to M* > 10 10 M Q . The additional incompleteness due to our 
velocity dispersion criteria is less than 10 % for all M* > 10 10,3 M Q and n > 3, but is significant for n < 1 at all masses. 

There is one complicating factor in the comparison between this general field sample and the G09 sample that we discuss in 
the main text. For the field sample, we are forced to rely on either the De Vaucouleurs/exponential model fits provided in the 
basic SDSS catalog, or the Sersic fits given in the NYU VAGC. (Recall that we need a measure of total magnitude to derive 
M*, and both an effective radius and a Sersic index measurement to derive Md,„.) Both of these sets of measurements have their 
faults. The SDSS model fits are overly simplistic in that they assume that n is equal to either 1 or 4; this will clearly introduce 
systematic errors in the fit quantities a s a function of (intrinsic) profile shape. The VAGC Sersic fits are also known to suffer 
systematic errors (Blanton et al. 2005a), due to background over-subtraction (G09). 

Bearing both these issues in mind, in Figure[TT]we show the relation between M„ and „ for our general, field galaxy sample, 
using either the Sersic fit parameters from the NYU VAGC (left panel) or the model fit parameters from the basic SDSS catalog 
(right panel) . Using either set of parameters, the results for this general sample agree really very well with what we have found 
for the G09 sample in Figure |4] Moreover, the two panels in Figure QT| agree remarkably well with one another, even despite the 
significant and very different systematic errors that each set of measurements suffers from. 

How can this be? It turns out that the covariance between Sersic parameters leaves ratio M»/Md.„ is remarkably robust to 
both random and systematic errors in the Sersic fits, provided „ is calculated as per Equation [2] To illustrate this, let us 
compare the IG09I and VAGC measurements. Although there are large differences in all three parameters individually, there 
are tight c orrelations betwee n An, Am tot , and AR e . (Here and in what follows the 'A' implies the difference between the 
VAGC and lGou et al.l (l2009l) -derived value, in the sense of VAGC-minus-G09.) Now, M* scales directly with total flux; fitting 
to AM* as a function of An, we find that AlogM* oc 0.04An, with an rms scatter in AlogM* of 0.07 dex. M^„, at least as 
defined in Equation [2] depends on both the effective radius and Sersic index. Considering the change in „ due to changes 
in size alone, we find AlogMd,„ oc 0.11 An, with a scatter of 0.07 dex; for the effect due to changes in the Sersic index alone, 
we find AlogMd,,, oc -0.06 An, with an rms scatter of 0.03 dex. Taken together, the overall change in Md.„ scales with An as 
AlogM d ,„ oc0.04An. 

Thus we see that the changes in M* and M^„ thus have virtually the same dependence on An, leaving the ratio M*/Md.„ 
virtually unchanged. Further, the scatter in A(M»/Md.„) is just 0.04 dex. Using the basic SDSS model fits, while we find slightly 
stronger dependences with An, we still find that the ratio M*/Md,„ remains very robust. We stress that the above argument only 
holds if we account for the dynamical effects of structure in the calculation of M±„: if we removed the n-dependence of Md,„ 
that enters via Kv(n), then we would find that the ratio Alog(M*/Md) oc -0.06An, in agreement with the expectation from the 
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FIG. 12. — Comparing dynamical and stellar mass for the G09 sample, using structural fit parameters from the NYU VAGC or the basic SDSS catalog. — All 
symbols and their meanings are directly analogous to Figures Bland II II In comparison to Figure RJ this Figure differs only in that we have used the Sersic-fit 
parameters from either the NYU VAGC (Blanton et at. 2005a) or the De Vaucouleurs/exponential model fit parameters from the basic SDSS catalog. The fact 
that these results agre e ve ry well with those shown in FigureRjshows that the ratio M* /M&„ is not extremely sensitive to errors in the structural fit parameters. In 
comparis on t o FigurefTT] the difference is that we only show galaxies in the G09 sample. The fact that the results in this Figure agree very well with those shown 
in Figure [TTI show that sample selection effects do not play an important role in our results. 




FIG. 13. — The M*-Md „ relation for galaxies in different states of activity. — Each panel of this Figure shows the relation between M» and „ for different 
subsamples of the general 0.035 <z< 0.08 galaxy population, split according to their spectral classification, and analyzed using the Sersic fits given in the NYU 
VAGC. From left t o righ t, we show non-emission galaxies, star forming galaxies, and AGN/composite spectra galaxies; the spectral classifications are those of 
Brinchmann et al. (2004), which are based on the BPT diagram. In each panel, the heavy dashed line shows the best fit M*-Md „ relation; for comparison, the 
grey short-dashed line shows the best fit relation for the sample as a whole. While each subsample follows a similar M*-Mj „ relation, there is an offset between 
the different relations, on the order of 0.07 dex. As for the G09 sample, the differences between the different subsamples disappears if we consider only n > 3 
galaxies. That is, these offsets appear to be due to the different distributions of n within each subsample, rather than intrinsic differences in the values of the 
stellar-to-dynamical mass ratio for galaxies in different states of activity. This argues the idea that selection effects play a major role in shaping our results. 



analysis immediately above. 

To explicitly demonstrate that the observed relation between M* and M±„ is not particularly sensitive to the measurements used 
to derive the values of M* and M±„, in Figure[l2]we show the M*-M±„ relation for the IG09I sample analyzed using structural 
parameters from the NYU V AGC (left panel) or from the SDSS catalog (right panel). In comparison to Figure |4] the slo pe of 
the M*-M d „ relation for the|G09 sample is very similar using any of the three sets of structural parameters: 0.92 for the G09 
fits, compared to 0.89 for the NYU VAGC fits, and 0.88 for the SDSS model fits. The normalization of the M*-M dj „ relation is 
slightly more sensitive: at M di „ = 10 11 M Q , we find that Alog(M*/M djn ) = -0.23, -0.24, and -0.32 dex using the lG09i VAGC, 
and SDSS fits, respectively. 

In comparison to Figure [TT] the results in Figure Q~2] also demonstrate that the M»-M d „ relation for the IG09I is very similar to 
that for a more general field galaxy sample. For example, using structural parameters from the VAGC, the logarithmic slope and 
intercept of the M»-M d „ relation are a = 0.89 and bn = -0.24 for the 1G09I sample, compared to a = 0.91 and b\i = -0.27 for the 
general galaxy sample. 

In Figure [13] we show that galaxies in different states of activity follow very similar M*-M d „ relations. In this Figure, we 
have split the general ga laxy sample into non-emission, star forming, and AG N and composite spectra subsamples using the 
Brinchma nn et al.l (12004) spectral classification scheme described in Section l5~2l . The logarithmic slopes of the M»-M d „ relation 
for each subsample agree with one another, as well as with that for the sample as a whole, to within a few percent. We do find that 
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FIG. 14. — Separating out the mass- and structure-dependence of the mass ratio M*/Mj „ of a field sample of 0.035 < z < 0.08 galaxies, showing active and 
passive galaxies separately. — In analogy to Figure|6] the colored lines in the left panels show the mass-dependence of the mass ratio M»/Md,„ in bins of Sersic 
index; those in the right panel shows how M* /M$ „ varies with n in bins of M* . The precise bins are given within each panel, along with the best fit parameters for 
the relation for each bin. The results in this panel have been derived using the NYU VAGC Sersic fit parameters. In comparison to Figure[6] we see qualitatively 
similar behavior for M* /M±„ as a function of both M, and n for the general galaxy sample as we do for the heavily-selected G09 sample. Moreover, we point out 
that this is true for non-emission, star forming, AGN and composite spectra galaxies separately, as well as for the general sample as a whole. Because the NYU 
fits suffer systematic biases, there are quantitative differences in the results shown in this Figure and those in Figure [6] Despite these quantitative differences, 
the relatively weak dependence of M* /M±„ on n for each bin in M, supports our main result; viz., that accounting for non-homology leads to reasonably good 
consistency between stellar and dynamical mass estimates. 



the M»-Md.„ relation for the star forming subsample is offset from th at fo r the non-emission and AGN/composite subsamples, at 
the level of 0.07 dex. However, as for the IG09I sample (see Section IBT2I 1. these small differences disappear if we consider only 
n > 2.5 galaxies. We thus conclude that these differences are principally driven by the different distribution of n values within 
the star forming sample, rather than intrinsic differences in the stellar-to-dynamical mass ratios of star forming galaxies. 

Finally, in FigureO we separate out the M* - and n-dependences of M* / Md,„ for the general galaxy sample, analyzed using the 
NYU VAGC Sersic-fit parameters. Again, we find that the ratio M*/M± n depends on both (at fixed n) and on n (at fixed M*). 
The results in this Figure suggest that the mass-dependence of M*/M±„ may flatten considerably for n < 2 and 10 < logM*/M Q 
< 10.5 (ie. below the mass limit of the G09 sample). 

As we have noted above, there are significant differences in the values of M* and n given by G09 and in the VAGC. These 
differences in M* and n mean that the results in this Figure are not in quantitative agreement with those shown in FigureQ4] even 
though both datasets show good agreement in the global M*-Md,„ relation. In particular, the slope of the M*/M± n -n relation at 
fixed is significantly steeper than we find for t helG09l sample. This is at least partially due to the bias in the NYU values of n; 
the NYU values are systematically lower than the G09 values, which has the effect of steepening the M*/M± n -n relation. 

With these caveats, the main conclusion to be drawn from Figure [14] is that accounting for non-homology in the derivation 
of dynamical masses leads to considerably better consistency between and M± n (as a function of «, and at fixed M*), in 
agreement with our findings in Section @] Further, we note that we find similar and consistent behavior in M*/Md,„ at fixed M* 
and n for each of the three subsamples shown in FigureQj] in agreement with our conclusions above. 

In summary, then, in this Appendix we have demonstrated two things. First , we have shown that we find very similar results 
for the G09 sample, analyzed using the results of the Sersic fits given by G09, and for a more general galaxy sample, analyzed 
using either the Sersic fits given in the NYU VAGC or the De Vaucouleurs/exponential model fits given in the basic SDSS 
catalog. Secondly, we have shown that we find very similar results for the G09 sample analyzed using any of these three sets of 
structural parameters. The most important conclusion to be drawn from these results is that the results we have presented in the 
main text are not driven, nor particularly sensitive to, selection effects. 
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